Transition from diffusive to ballistic dynamics for a class of finite quantum models 
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The transport of excitation probabilities amongst weakly coupled subunits is investigated for a 
class of finite quantum systems. It is demonstrated that the dynamical behavior of the transported 
quantity depends on the considered length scale, e. g., the introduced distinction between diffusive 
and ballistic transport appears to be a scale-dependent concept, especially since a transition from 
diffusive to ballistic behavior is found in the limit of small as well as in the limit of large length 
scales. All these results are derived by an application of the time-convolutionless projection operator 
technique and are verified by the numerical solution of the full time-dependent Schrodinger equation 
which is obtained by exact diagonalization for a range of model parameters. 

PACS numbers: 05.60.Gg, 05.70.Ln, 05.30.-d 



o 

(D 

e 

I 



O 
o 



> 

O 
p 

o\ 
o 
1^ 
o 



X 



There are certainly many ways in order to analyze the 
transport behavior of quantum systems: Linear response 
theory such as implemented in the Kubo-formula [l|, 0, 
Si Si Si @| 7 mcthods which intend to map the quantum 
system onto a system of interacting, classical, gas-like 
(quasi-)particles 7 . iS. iol. [l^. [Til (possibly represented by 
Green's functions [ij, [iJDi methods which include 
reservoirs of the transgorted quantity, either explicitly 
[lil or effectively [16, 17, 18 1, and probably many more. 
However, except for those involving reservoirs, none of 
these methods allows to investigate the dependence of 
transport behavior on the length scale of a finite system. 
Unfortunately, the remaining methods which do involve 
baths are conceptually more subtle and computationally 
rather challenging Tq, 17, 18| . 



In this letter we will discuss a class of finite models 
which indeed exhibits scale-dependent transport of, e. g., 
excitation probabilities, energy or particles, according to 
the respective interpretations of the model. In particular, 
we will develop a method which allows for the complete 
characterization of the available types of transport and 
their dependence on the considered length scale. It will 
be especially demonstrated that the numerically exact 
solution of the full time-dependent Schrodinger equation 
perfectly agrees with the theoretical predictions which 
are obtained from this method: For intermediate length 
scales diffusive behavior is found, that is, the dynamics of 
the transported quantity is well described by a respective 
diffusion equation and the spatial variance of an initial 
excitation profile increases linearly in time. Contrary, in 
the limit of small as well as large length scales diffusive 
behavior breaks down. This non-diffusive transport type 
turns out to be "ballistic" which is to be understood in 
terms of a quadratic increase of the spatial variance. (See 
the last page of this letter for a classification of these 
findings in the framework of standard solid state theory.) 

According to Fig. [l] our model consists of N identi- 
cal subunits which are labelled by fi — 0,1, . . . , N — 1. 
These subunits are assumed to have a non-degenerate 
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FIG. 1: A model of A'^ identical, weakly coupled, subunits 
featuring a non-degenerate ground state, an energy gap AE 
and an energy band 5e with n equidistant states. The dots 
indicate excitation probabilities and are supposed to visualize 
a state from the investigated "single-excitation space" . 



ground state with energy eg = as well as n excited 
states with equidistant energies Si, i — 1,2, ... ,n. We 
will focus on the (invariant) "one-excitation subspace" 
which is spanned by the set where \fJ.,i) corre- 

sponds to the excitation "sitting" in the ith exited state 
of the /ith subunit. Using this notation, the total Hamil- 
tonian H = Hq + XV, the sum of a local part Hq and an 
interaction part V, is given by 

^0 = ^^<^i IMi«)(Mi«I I (1) 

V = ^^q,|mii)(m+1ii| + H.c. (2) 



We use periodic boundary conditions, i. e., we identify 
fi = N with fj, = 0. The Cij form a normalized Gaussian 
random matrix with zero mean, A is an overall coupling 
constant. Obviously, V corresponds to nearest-neighbor 
hopping. 

This system may be viewed as a simplified model for, 
e. g., a chain of coupled molecules or quantum dots, etc. 
In this case the hopping of an excitation from one subunit 
to another corresponds to transport of energy, especially 
if AE ^ Se. The model may also be viewed as a tight- 
binding model for particles on a lattice. In this case the 
hopping corresponds to transport of particles. There are 
n "orbitals" per lattice site but no particle-particle in- 
teraction in the sense of a Hubbard model [cf. Eq. 
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Due to the independence of Cij from /x, these are sys - 
tems without disorder in the sense of, say, Anderson [l9| , 
in spite of the fact that the Cy are random. For some 

Hterature on this model class, see [1,0, mm. 

The total state of the system is naturally represented 
by a time-dependent density matrix p{t). We denote by 
= l/i, i) (ij, i\ the operator which projects onto the 
exited states i) of the /xth subunit's band. Thus, the 
quantity Pfj,{t) — Tr[n^p(i)] represents the probabihty 
for finding an excitation of the /ith subunit, while all 
other subunits are in their ground state. 

Our aim is to describe the dynamical behavior of these 
local probabilities and to develop explicit criteria that 
enable a clear distinction between diffusive and ballistic 
transport. We will call the behavior diffusive, if the P)j,{t) 
fulfill a discrete diffusion equation 



Pn 



-2P„, 



(3) 



with some diffusion constant k. As well known from the 
work of Fourier, such a diffusion equation decouples with 
respect to, e. g., cosine-shaped spatial density profiles, 
i. e., the above equation yields 



Fq^-2{l-cosq)KFg , Fy 



CgY.C08{qfl)P^ (4) 



with q = 2tt k/ N, fc = 0, 1, . . . , N/ 2 and an appropriate 
normalization constant Cq. Hence, if the model indeed 
shows diffusive behavior, then all modes Fq have to relax 
exponentially. But if the modes Fq are found to relax 
exponentially only for some regime of q, the model is said 
to behave diffusively on the corresponding length scale 
I = 2tt/ q. The dynamics of the Fq, as resulting from 
the quantum system, is most conveniently expressed in 
terms of expectation values of "mode operators" 



F,(i) = Tr[$,p(i)] , ^q = CqY,C08iqfi)n^, (5) 



where we now choose Cq = ^j2lnN for g 7^ 0, tt as well 
as Co = = yJl/nN such that Tr[ = 5qq>, the 
mode operators are orthonornial. 

A strategy for the analysis of the dynamical behavior of 
the Fq is provided by the projection operator techniques 
0, 13 US 0| ■ III order to apply these techniques one 
first defines an appropriate projection superoperator V . 
Formally, this is a linear map which projects any density 
matrix p{t) to a matrix V p{t) that is determined by a 
certain set of relevant variables. Moreover, the map has 
to be a projection in the sense oiV'^ = V. In the present 
case a suitable projection superoperator is defined by 

V p{t) ^ J2 % Pit)] ■ (6) 

Due to the orthonormality of the mode operators this 
map is indeed a projection operator. The great advantage 



of this approach is given by the fact that it directly yields 
an equation of motion for the relevant variables, in our 
case the Fourier amplitudes Fq. We concentrate here on 
a special variant of these techniques which is known as 
time-convolutionless (TCL) projection operator method 



2a,l27[. Considering initial states with V p{0) — p(0), 



this method leads to a time-local differential equation of 
the form Fq{t) — ~rq{t) Fq{t). So far, this formulation 
is exact. Modes with different q do not couple due to the 
translational invariance of the model. 

Nevertheless, the exact Fq{t) is hard to determine but 
the TCL method also yields a systematic perturbation 
expansion for the rate ^q(t) in powers of the coupling con- 



stant A, namely, Tqit) = X^rq^2it) + X'^Ti^qit) 



(all 



odd order contributions vanish). Note that a truncation 
of the above TCL expansion may yield reasonable results 
even and especially in the regime where the elements of 
the perturbation matrix are larger than the level spacing 
oi Hq, i. e., the regime where standard time-independent 
perturbation theory breaks down. 

In leading order a straightforward calculation yields 



Fq{t) = -2{l-cosq)-f2it)Fq{t) 



(7) 



where the rate 72 (i) is defined by 72(0 = /g dr fir) and 
the two-point correlation function /(r) is given by (time 
arguments refer to the interaction picture) 



/(r) = -Tr[T/(t)T/(ti)n, 



t - U 



This function is completely independent from p, because 
of the translational invariance of the model. Of course, 
/(r) depends on the concrete realization of the random 
interaction V . But due to the law of large numbers, the 
crucial features of V are nevertheless the same for the 
overwhelming majority of all realizations of V , as long as 
y/n 3> 1. And in fact, /(r) typically assumes the form 
in Fig. O It decays like a standard correlation function 
on a time scale of the order Tc — 1/ 5e. The area under 
this first peak is approximately given by 7 = 27r n A^/ 5e. 
However, unlike standard correlation functions, due to 
the equidistant level spacing of the local bands, /(r) is a 
strictly periodic function with the period T = 27rn/ 5e. 
Consequently, its time-integral 72 (i) nearly represents a 
step function, see Fig. [2] 

Thus, for Tc < t < T we find from Eq. ^ 



Fq{t) = -2 {1- cos q)^Fq{t) 



(8) 



The comparison with Eq. ^ clearly indicates diffusive 
behavior with a diffusion constant /t = 7. And indeed, 
for modes which decay on a time scale t with Tc < t < T 
we find an excellent agreement between Eq. ([8]) and the 
numerical solution of the full time-dependent Schrodinger 
equation, obtained by incorporating Bloch's theorem and 
exactly diagonalizing the Hamiltonian within decoupled 
subspaces. A "typical" example for a single realization of 
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FIG. 2; Sketch of the correlation function /(r) [dashed line] 
and its integral 72 (t) [continuous line]: /(r) features complete 
revivals at multiples of T such that 72 (t) has a step-like form. 
This sketch indeed reflects the numerical results for /(r) and 
72(4) but highlights the relevant time scales for clearness. 



the random matrix V is shown in Fig. [3] Note that the 
size of the system is chosen adequately large such that 
there are only small fluctuations of the order of 1/n which 
arise due to the discrete spectrum, see, e. g., 20l. 21 1. 




time t [T] 



0.08 



FIG. 3: Evolution of a Fourier mode Fjyit) which decays on an 
intermediate time scale Tc < t <T. Numerics (crosses) shows 
an exponential decay which indicates diffusive behavior and 
is in accord with the theoretical prediction (continuous line). 
Parameters: = 120, n = 500, Se = 0.5, A = 0.0005. 

However, until now the above picture is not complete 
for two reasons. The first reason is that Fq{t) may decay 
on a time scale that is long compared to T. According 
to Eq. Q, this will happen, if 2(1 - cosq)^T > 1 is 
violated. If we approximate 2(1 — cosq) k, = 4tt^/P 
for rather small q (large /), this leads to the condition 



I5e 



> 1 



(9) 



If this condition is satisfied for the largest possible I, i. e., 
for ^ = A^, the system exhibits diffusive behavior for all 
modes. If, however, the system is large enough to allow 
for some I that violates condition ([5]), diffusive behavior 
breaks down in the long- wavelength limit. This result is 
again backed up by numerics, see Fig. [4l 

Towards what transport type does the system deviate 
from diffusive, if condition ((9]) is violated? In this regime 
we have to consider time scales with t 3> T, as already 
mentioned above. We may thus approximate 72(0 ~ 
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FIG. 4: Deviations of the time evolution of the Fourier modes 
Fq{t) with q = 2n/N, the longest wavelength, from a purely 
exponential decay for different model parameters A'' and n. 
These deviations are based on a measure used in [22] and are 
in accord with the claim that diffusive transport behavior is 
restricted to the regime defined by condition ([9]). Other model 
parameters: Se = 0.5, A = 0.0005. 



27t/T, see Fig.H Plugging k ^ 2jt/T into Eq. ^ yields 
a spatial variance ct^ = (^^) — (/i)^ = 2-ft'^/T, contrary 
to = 2^t which results for k = 7. (k = 7 applies in 
the regime where condition ([9]) is satisfied.) This change 
of the time-scaling of clearly indicates the transition 
from diffusive to ballistic transport. The validity of the 
TCL approach in the ballistic regime is again backed up 
by numerics: Here, the TCL theory obviously predicts a 
Gaussian decay of Fq{t) which is in agreement with the 
numerical results, see Fig.[S] 
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FIG. 5: Evolution of a Fourier mode -F,r(i) which decays on 
a time scale t ^ Tc, parameters: = 120, n = 500, 5e = 0.5, 
A = 0.004. Inset: Evolution of a Fourier mode -F,r/6o(^) which 
decays on a time scale t^T, parameters: A'^ = 120, n — 500, 
Se = 0.5, A — 0.0005. In both cases numerics (crosses) shows 
a Gaussian decay which indicates ballistic behavior and is in 
accord with the theoretical predictions (continuous lines). 

In a second case transport may be non-diffusive, if the 
corresponding Fq(t) decay on a time scale that is short 
compared to Tc- This will happen, if 2 (1 — cos q)^Tc^ 1 
is violated. We may approximate 2(1 — cos g) « 4 for the 
largest possible q (smallest possible I). Hence, the above 
inequality may be written as 



(5e2 



< 1 



(10) 



If this inequality is violated, diffusive behavior breaks 
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down in the limit of short-wavelength modes. Moreover, 
if the second order still yields reasonable results for not 
too large A, we expect a linearly increasing rate 72 {t) and 
thus a Gaussian decay, that is, according to the above 
reasoning, ballistic transport. For increasing wavelength, 
however, the corresponding inequality will eventually be 
satisfied, thus allowing for diffusive behavior. Also these 
conclusions are in accord with numerics, see Fig. [5l 

Standard solid state theory always predicts ballistic 
transport for a translational invariant model without 
particle-particle interactions. Nevertheless, in the limit 
of many bands (many orbitals per site) and few sites (few 
fc-values) two features may occur: i) The band structure 
in fc-space becomes a disconnected set of points rather 
than the usual set of distinct smooth lines. It is hence 
impossible to extract velocities by taking the derivatives 
of the dispersion relations, ii) The eigenstates of the 
current operator no longer coincide with the Bloch eigen- 
states of the Hamiltonian and the current becomes a non- 
conserved quantity, even without impurity scattering. It 
is straightforward to check that both features occur in 
the regime where condition ^ is fulfilled. This is the 
regime where standard solid state theory breaks down 
due to the fact that the system is to "small" . 

On the other hand, our numerical simulations clearly 
reveal that the TCL projection operator technique is ap- 
plicable both in the diffusive and in the ballistic regime 
and that it correctly describes the transition between 
these regimes. Recall that all analytical results have been 
obtained from the second order contribution [Eq. ([7])] 
of the TCL expansion. And in fact, it is possible to 
demonstrate that higher-order contributions are confi- 
dently negligible. We omit the details of the determina- 
tion of the higher orders of the TCL expansion (several 
examples are discussed in (26l|). since this is beyond the 
scope of this letter and since the excellent agreement with 
the numerical simulations is surely evident. 

However, the fact that already the second-order TCL 
equation ([7]) yields an excellent quantitative agreement 
with the numerics, by contrast to standard perturbation 
theory, can be understood from the following argument. 
We first note that for A = the energy eigenvalues of the 
Hamiltonian are A^-fold degenerate, corresponding to the 
resonant transitions — > |/i± between neighbor- 
ing subunits induced by the interaction. While standard 
perturbation theory is spoiled by the presence of these 
resonant transitions, the TCL expansion works perfectly 
well even in the case of exact degeneracies. In fact, an 
analysis of the structure of the higher-order correlation 
functions reveals that for an interaction Hamiltonian Vd 
which includes only the resonant transitions, the second- 
order TCL equation ([7]) is exact on average (for arbitrary 
A) with small fluctuations of order l/n. This is due to 
the fact that the TCL method is based on an expan- 
sion in terms of the ordered cumulants of Vd [2^, the 
higher-orders of which vanish by virtue of the Gaussian 



character of the coupling matrix elements. This implies 
that the second order TCL equation already reproduces 
exactly the effect of all resonant transitions, which is the 
crucial reason for the success of the TCL projection op- 
erator approach. 
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